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^ ■ A result of Zyczkowski and Sommcrs [J. Phys. A 33, 2045-2057 (2000)] gives the eigenvalue 

Ph. probability density function for the top N x N sub-block of a Haar distributed matrix from 

. U{N + n). In the case n > N, we rederive this result, starting from knowledge of the 

distribution of the sub-blocks, introducing the Schur decomposition, and integrating over all 
variables except the eigenvalues. The integration is done by identifying a recursive structure 
which reduces the dimension. This approach is inspired by an analogous approach which has 
been recently applied to determine the eigenvalue probability density function for random 
matrices A^^B, where A and B are random matrices with entries standard complex normals. 
We relate the eigenvalue distribution of the sub-blocks to a many body quantum state, and 
to the one-component plasma, on the pseudosphere. 



■ 1 Introduction 

> 

■ The plane, sphere and pseudosphere are special geometries in a number of related many body statistical 
^SJ . systems. In particular the lowest Landau level wave function for quantum particles in a magnetic field, 

\ the equilibrium statistical mechanics of a one-component plasma at a special value of the coupling, 

and the eigenvalue probability density function for certain random matrix ensembles all give rise to a 
' special solvable state in these geometries. In regard to the first of these, regarding the geometries as 

■ surfaces in appropriate 3-spaces, and imposing a perpendicular magnetic field, the solvable state is the 
ground state wave function, in that it can be written in an explicit factorized (Jastrow) form and the 

I correlations can be expressed explicitly as determinants [H [5] . The absolute value squared of these wave 

^ ' functions have the interpretation as the Boltzmann factor for a one-component plasma, confined to the 

corresponding surface, and interacting at the special value of the coupling F q^/ksT = 2 [HITS]. And 
the Boltzmann factors, after projecting onto the plane via an appropriate stereographic projection in the 
case of the sphere and pseudosphere, allow for realizations as the eigenvalue probability density function 
for three particular random matrix ensembles |121 123|, I16j . The determinant form for the correlations 
in the quantum problem carries over to the plasma and eigenvalue distributions, so these many body 
systems are also solvable. Although not to be addressed here, we mention too that these same geometries 
play a special role in relation to the study of the zeros of random polynomials [HI HH] • 

It is the purpose of the present paper to contribute to the theory of the random matrix ensemble 
relating to the pseudosphere, or more precisely to the Poincare disk. Specifically, our interest is in the 
eigenvalue distribution of the ensemble of the top N x N sub-block of Haar distributed matrices in 
U{N + n) (due to the invariance of Haar distributed unitary matrices with respect to left and right 
multiplication by permutation matrices, in fact sub-blocks formed by the entries of any TV rows and N 
columns will have the same distributions as the top N x N sub-block). This ensemble has shown itself 
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in a number of recent works in random matrix theory [51 llll 1171 121] . We present a new derivation of 
a result due to Zyczkowski and Sommers [23], giving that the eigenvalue probability density function is 
proportional to 

N 

Y[{i-\zifr-'x\.A<i n i^fc-^.f' (1-1) 

1=1 l<j<k<N 

(here XA = 1 if A is true, and XA — otherwise). Our proof, although applying only to the case when 
n > TV, is arguably simpler than the original. Moreover, our derivation unifies the case of the sphere 
and pseudosphere, by showing that a recently formulated approach to the eigenvalue probability density 
function of the random matrix ensemble corresponding to the sphere [13] , can be adapted to the case of 
the pseudosphere. 

2 The pseudosphere and plasma system 

The pseudosphere refers to the two-dimensional hyperbolic space with constant negative Gaussian curva- 
ture K = — [1 . It is naturally embedded in the three dimensional Minkowski space with coordinates 
{yo,y 1,1/2) and line element 

(dsf = -{dyof + {dy,f + {dy^f. (2.1) 

The pseudosphere is defined as the upper branch of the equation 

2,2,2 2 
-2/0 + 2/1 + 2/2 -a 

and can be parameterized by 

yo = acoshr, = asinhrcos(/), j/2 = asinhrsin^. 

In terms of this parametrization the line element (j2.ip reads 

{dsf = a^{dTf + c? sinh^ r {d<\)f . (2.2) 

Furthermore, the geodesic distance between (r, </>) and (r', </)') is such that 

cosh — — cosh T cosh r' — sinhrsinhr' cos(0 — 0'), (2.3) 
a 

while the volume element corresponding to the surface area is 

dS ^ (? sa^T dTd4>. (2.4) 

With z := x+iy, the pseudosphere is projected onto the Poincare disk via the stereographic projection 

^ = 2atanh^e*'^, \z\ < 2a (2.5) 
(equivalently (2/0,2/1,2/2) ^ 2a (yi/j/o, 2/2/2/0)) and then 



dxdy 

(1 -'WJ^a^y 



ds= ^. (2.6) 
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Note that with z expressed in polar form, (12. 5p is equivalent to writing 



r 

2atanh — , \z\ < 2a. 



In terms of the coordinates (r, cf)) the Laplacian A is computed from (|2.2p as 

1 / 1 9 . , 9 1 
A — —{ — - — -— smh T- 



A 



a? \ sinh r dr dr sinh^ t i 

Alternatively, in terms of the coordinates (r, (p) the Laplacian reads 

From the Laplacian we can specify the Coulomb potential on the pseudosphere. In fact, two potentials 
$ and $ will be defined. With sjk denoting the geodesic separation between particles j and k, and A^ 
the Laplacian with respect to the particle j, the first is specified as the solution of the Poisson equation 

A,$(s,fe) = -2TrS{sjk) (2.8) 

subject to the boundary condition that <I>(s) ^ as s ^ oo, and the second as the solution of the Poisson 
equation 

A,$(s,fe) - -2nS{s,k) - (2.9) 

In both (|2.8p and (|2.9p 6{sjk) denotes the delta function relative to the volume form (|2.4p . The Poisson 
equation (|2.8p corresponds to the potential at {tj, due to a unit charge at (t^, (j>k), while (|2.9p corre- 
sponds to the potential at (t, , 0j) due to a unit charge at (r^, (/)fc) and a uniform smeared out charge of 
the same sign and density l/47ra^. The solution of (|2.8p is readily verified to be 



<!>(.,,) = -logtanhf . - log f^-^ , (2.10) 
and that of ([^ to be 

*(.,) = - ios...„hif = - ,og ( (i_(,,^,„^^^;'{';^^,„,.,v, )- 

Choosing one or other of (|2.10p . (|2.1ip . we want to form a one-component plasma system, consisting 
of N mobile unit charges, and a smeared out uniform background charge density of net charge —77. In 
the large N limit, we know from [15] that the bulk equilibrium statistical mechanics is the same for 
both. However, our interest is in finite TV, and to obtain an analogy with a random matrix eigenvalue 
probability density function it is necessary to base the plasma system on the potential (|2.1ip . Since 
(|2.1ip corresponds to (|2.9p . to get a net background charge density —t] we must also cancel the smeared 
out positive charge implied by the constant term therein, and thus impose a smeared out charge density 
—rj — N/Atto?'. This will interact with the particles to create a potential V{r), which must satisfy the 
Poisson equation 

AV{r) = 27r(7y + iV/47ra2). 

Choosing for convenience ^(0) — (the constant terms in the total potential are of no interest to us), 
and using the form (|2.8p of A, the solution of this equation is seen to be 

2 

V{t) = -(27177^2 ^ j^i^-^ (1 - ^) • (2-12) 



{zj - Zfc)/2a 
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The total potential energy U, excluding constants, is obtained by summing the potential (12. lip over all 



pairs, and adding to this ^('^i)- We thus have [15] 

N 

J2 log {\zk - z,\/2a) - (27rr;a2 + 1/2) ^log(l - r]/Aa^) 



U = - 

l<j<k<N j=l 



implying that the Boltzmann factor is such that 

N 

e-^^dSi ■■■dSNOZ Y[{1 - |^;|2/4a2)(2^r,a^ + l/2)/3-2 "Q _ f^Xidyi ■ ■ ■ dxNdyN{2.13) 

1=1 l<j<k<N 

3 Magnetic analogy 



It is remarked in [TJ] that with / :— y/hc/eB the magnetic length 

V^,(.)=(l-U_) J =0,1,... (3.1) 

are, up to normalization, the lowest Landau level eigenstates for a quantum particle on the pseudosphere, 
in the presence of a constant perpendicular magnetic field B, and projected onto the Poincare disk. 
Forming a Slater determinant out of the first N of these gives an A^-body state ip such that jV'l^ is 
proportional to (|2.13p with 

/? = 2, AiTTja^ + 1 = 2(a/0^ (3.2) 

We take this opportunity to gives some details on how (|3.ip comes about (see also [?])• 

Let r denote the unit vector perpendicular to the pseudosphere in Minkowski space. By the formal 
analogy with the surface of a sphere with radius R and polar co-ordinates (0, 4>) as specified by 

a <-» iR, T ^ iO 

one sees that the vector potential 

1= f 0,0,^^(1 -cosh r)) (3.3) 
V smh r / 

is such that V x A = Br and so corresponds to the sought perpendicular magnetic field. Further, it 
follows similarly that the Hamiltonian 



Im \ c / 



2 

2m 



is such that 



/ 1 (9 / . , d \ , 1 / d /a\2, 

H 



2ma^ V sinhrc^TV drJ sinh 

(of course in writing down p.3p and p.4p knowledge of the corresponding results for the sphere are 
assumed [22]). The change of variables (|2.5p so as to project to the Poincare disk gives that p.4p reads 

;j2 . , U|2, /^2 , U|2, , /5 a l^|2s 



2m 
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Following [5], with lOc eB /mc, we note the projected Hamiltonian (|3.5p is of the form 
with 



But (13.61) is such that 



and so we are lead to the result that any function of the form e~'^f{z) with / analytic in z belongs to 
the lowest Landau level (i.e. is an eigenfunction of H with eigenvalue hujc/i). To further specify /, one 
notes that 

d d 
oz oz 

commutes with H as defined by (|3.6p for general $ = $(|zp), and seeks simultaneous eigenfunctions of 
H and J. Such eigenfunctions are given by 

e-*z^ J =0,1,... 

and p.ip follows. 

Now if all the particles are fermions, the ground state ip is formed out of a Slater determinant of the 
single particles states, 

. . . , Zat) = ^det[V'j-l(Zj)]j,fc=i,...,Ar 



^n(i-^)'°"'" n ft-'.) (") 



/=1 i<j<k<N 

where C is an appropriate normalization, and the second equality follows upon use of the Vandcrmonde 
determinant identity. The modulus squared is indeed identified with the case (I3.2p of (|2.13p . 



4 Random matrix ensembles 
4.1 The spherical ensemble 

The eigenvalue probability density function for matrices Y = A^^B with A and B random N x N 
matrices consisting of independent complex Gaussian entries, turns out to be a one-component plasma 
on the sphere. This spherical analogue of (jl.ip was proved in |17| . Below, we give a simplified and unified 
approach that works for the sphere as well as the pseudosphere. We first sketch the proof for the sphere 
(more details may be found in the forthcoming book [13j). The first step is to determine the distribution 

of r. 

Proposition 1. The distribution ofY is proportional to 



det(iAr + yyt)2w 
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Proof. With X a complex matrix and {dX) denoting the wedge product of the differentials of the real 
and imaginary parts of its elements, by noting 

{dB) = \detA\^^{dY) (4.2) 

it follows that the joint density of A and Y is proportional to 

|det^|2^e-^^^^P"+^^')-4'>. 

The distribution of Y is obtained by integrating this over A. For this we change variables C — A{1m + 
yyt)i/2 g^j^(j i^y g^j^ approprlatc analogue of (|4.2[) and integration over C, obtain the density of Y to be 

det(C)p^e-T^(gg') dC 
det (Ijv + y r t ) A' det (Ijv + YY'f ) ^ ' 

from which (|4.ip results. □ 

The second and final step is to use knowledge of the distribution (|4.ip to proceed to compute the 
eigenvalue probability density function. 

Proposition 2. jlb^ The eigenvalue probability density function for Y is proportional to 

^ 1 

IlTTTm^ n (4-3) 

;=i ^ l^'l > i<j<k<N 
Proof. According to the Schur decomposition, Y can be written 

Y ^ UTnU~^ (4.4) 

where [/ is a unitary matrix, and T^r an upper triangular matrix in which the diagonal entries, denoted 
{zj}j=i,...,Af, are the eigenvalues of Y . Furthermore, we know that (see appendix 35 of |19j or section 6.3 
of [13]) 

{dY)= [| \zk- Zj^^iU^dU) hf^^dzj{dfN) (4.5) 

l<j<k<N 

where TV denotes the strictly upper portion of Tjv. Noting from (|4.4p that YY'^ = UT^Tj^U'^ shows that 
(|4.ip is independent of U . The factorization of the U dependence in (|4.5p then implies that the eigenvalue 
distribution is obtained from (|4.ip by integrating out the dependence of TV and is thus proportional to 



n |2fe -Zjf TAr(zi,...,ZAr), I n{zi, . . . , zn) ■= j —— ^— j— (dfv). (4.6) 

Following [13] ■ we will compute the integral In by deducing a recurrence for the integrals 

T„ p(zi, . . . , z„) := / r— {dfn) (4.7) 

J det(I„ +T„T^)P 

where p > n (this ensures convergence). Analogous strategies can be found in [TU [SI [TU]. For this 
purpose, let u denote the last column of T„. We can then write 



-1 + Tn^iTl_^ + UU^ 



ZnU 



(4.8) 
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showing 

det(I„ + T„Tl) = (1 + |z„n dot (l„_i + Tn-iTl_^ + ) . 

V J- r I I ^ 

Noting that the final term in the determinant is a matrix of rank 1, this latter formula can be further 
reduced to read 

det(I„ +r„T^) - (1 + |z„ndet(I„_i +T„_ir,l_i)(l + ] u\ln-i +T,,^iTl_^r^u) . (4.9) 
Hence, by first writing 

In,p{zi, . . . , Z„) = / (dTn-l) / {du) 



det(I„ + TnTly 



we can substitute (14.91) to obtain 



/„,p(zi, . . . , z„) = 'I / (dTn-l) 



(1 + |z„|2)p y det(I„_i + T„_ir,ti)P 

X / (du) (l + -^-l—u\l,,_, + Tn-iTl_,y'u) (4.10) 



Making the change of variables 

v = (1 + |z„n-i/2(i„_i + r„_ir,ti)-i/2, 

we see from this that 



In,pizi, ■ ■ ■ , Zn) — 4- |z |2^p-ri+l ^"-l'P-l(^l' ■ ■ ■ ' ^"-l) 



(4.11) 



where 

a. = / j^^. (4.12) 

7 (i + tnii)p 

Iterating (|4.1ip with rt = 2Af, p = N shows 



^ 1 

/Ar(zi, . . . , Zat) CX [| (1 + |Z,P)^+1 



and this substituted in (|4.6p gives (|4.3p . □ 

We remark that the name spherical ensemble comes about because upon making the stereographic 
projection from the plane to the sphere, 



(cf. (|2.5p ) transforms (|4.3p to the form 



z = 2i?tan-e*'^ 
2 



l<j<k<N 



where fj is the vector in M.^ corresponding to the point (i?, 0j , (j)j ) on the sphere. 
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4.2 Truncated unitary matrices — the pseudosphere ensemble 

Working now with truncated unitary matrices, we would like to repeat the steps which lead to establishing 
the eigenvalue probability density function (|4.3p . The first task then is to specify the distribution of the 
top N X N sub-block of Haar distributed matrices from U{N + n). Our working requires that n > N. 

Proposition 3. Let the top N x N sub-matrices in question be denoted Qn- The probability density 
of Qn is proportional to 

dct(lAr-QjvQjv)""'^, (4.13) 
supported on QnQ\ < 1 (i-S- matrices Qn with all singular values less than one). 

Proof. Following [8], we give a variation of the proof given in f^. The proof relies on a standard 
result in random matrix theory giving a construction of a particular Jacobi ensemble in terms of Wishart 
matrices [201 13 • With c and d being independent N x N and n x N complex Gaussian random matrices 
respectively (all elements i.i.d. standard complex normals), let C — c^c, D — d^d. We require the fact 
that for n > the probability density of 

J := {C + D)-^^^C{C + D)-^^^ (4.14) 

is proportional to 

det(Iw - J)""^, (4.15) 

supported on J < 1. 

Let U denote an {N + n) x [N + n) Haar distributed unitary matrix. Let W denote the N x {N + n) 
sub-matrix formed from the first N rows, and W denote the n x {N + n) sub-matrix formed from the 
final n rows, so U has the block decomposition 



U ■ 



w 
w 



Further, let X be an (A^ + n) x N Gaussian matrix of independent standard complex normals. From the 
facts that WX is then also distributed as an x Gaussian matrix of standard complex normals, and 
WX as an n X A^ Gaussian matrix of standard complex normals, and that WX, WX are independent, 
we can take c — WX and d — WX in ()4.14|) to conclude from the above result that 

J = (X^ X)-^/^X''W''WX{X'' X)-^/^ (4.16) 

has probability density (|4.15p . 

Introduce the singular value decomposition by writing 

X = U1AU2 

for Ui an (A^ + n) x [N + n) unitary matrix, U2 a,n N x N unitary matrix, and A an (A^ -I- n) x A" matrix 
with diagonal entries (that is entries {i,i), for i < N) equal to the positive square root of the eigenvalues 
of XtX. We then have 

XiX^X)-'/^ = [Ui]N+n,NU2 
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where the notation [Z]p,q denotes the top p x (7 sub-block of Z. Using this in (|4.16p shows 

N+n.,NU2- 

But the distribution of J is unchanged under conjugation by unitary matrices, and in particular U2- 
Applying this conjugation, we see that 

J has the same distribution as \u\]N,N+nW^W\Ui]N+n,N- (4-17) 

Since Ui and U are independent, and U has Haar distribution, UUi has the same distribution as U and 
hence, W has the same distribution as VFt/i_2], where we have written Ui,i — \U]N+n,N and 

Ui^2 denotes the last n columns of U . From this, it immediately follows that [VF^VFJtv.jv has the same 
distribution as \U\]N^Nj^nW^W[Ui]M+n,N ■ Together with (|4.17p . we conclude that [W^VF]jv,JV has the 
same distribution as J. 

This tells us that the probability density of Q\^Qn is proportional to (|4.13p . The probability density 
of Qn itself, being the top N x N block of [/, [/ G U{N + n), must be invariant under conjugation by 
S G U{N) and under complex conjugation, and so be a function of (the trace of) QI^Qn- But in general 
for a p X q {p > q) complex matrix Y with the probability density / of the form / = f {Q\^Q n) , the 
probability density for Y^Y is proportional to {det Y^ Y)p~'^ f{Y^Y). Here p = q = N, so we conclude 
that (|4.13p holds too for the probability density of Q^v- D 

Knowing that Qjv is distributed according to (|4.13|) . the task now is to introduce the Schur decom- 
position of Qat, and integrate out all variables but the eigenvalues. This can be done using the same 
recursive reduction of the dimension of the integral as used in the proof of Proposition [21 and the sought 
probability density function thus obtained. 

Proposition 4. Let n > N. Then, the eigenvalue probability density function of Qn is proportional to 

Proof. Because Qat is a sub-block of a unitary matrix, all eigenvalues have modulus less than one. 
Introducing the Schur decomposition (|4.4p of Qn in (|4.13p . and recalling the change of variables formula 
(j4.5p . we see the eigenvalue probability density function is proportional to 

II |zfe -z,f Jjv(zi,...,2iv), Jn {zi, zn) -.^ f dct{lN -TnTI)"-^ {dTN), (4.18) 

where again Tjv denotes the strictly upper triangular portion of Tjv. 
To compute Jn by recurrence we introduce, for p > 

J„,p(zi, . . . , z™) := / det(I„ - Tm,Tl)P {df^) (4.19) 



(cf. (|4.7p ). Analogous to (|4.8p . with u denoting the last column of T^r, we can write 



Working now as in deriving (|4.9p from (|4.8p shows 
dct(I^-r^T^) - (l-|z^|')det(I,^-i-T,^-irl_0(l- _] u\l,^.i-T,,,.^Tl_^)-^u). (4.20) 
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Substituting (|4.2Qp in (|4.19p and proceeding now as in the derivation of (14. lip from (|4.10p shows 

^m^pi^l 1 • ■ • 1 -^rn) ^m,p(f \ ) ^ Jm—l.p-}-! {'^li • ■ • ■^m — l') 

where Cm,p is again given by (|4.12p . Iterating this equation m times, we get 

m 

J,nA^U Zm) = - \zkfr+P-\ (4.21) 

fc=l 

Applying (|4.2ip with m — N,p — n — N shows 

N 

jjv(zi,...,z^)cxn(i-kii')""' 

1=1 

and this substituted in (|4.18|) gives (fTTj) . □ 
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